Dipole-Quadrupole dynamics during magnetic field reversals 
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The shape and the dynamics of reversals of the magnetic field in a turbulent dynamo experiment 
are investigated. We report the evolution of the dipolar and the quadrupolar parts of the magnetic 
field in the VKS experiment, and show that the experimental results are in good agreement with 
the predictions of a recent model of reversals: when the dipole reverses, part of the magnetic energy 
is transferred to the quadrupole, reversals begin with a slow decay of the dipole and are followed 
by a fast recovery, together with an overshoot of the dipole. Random reversals are observed at the 
borderline between stationary and oscillatory dynamos. 
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Despite a large variability of the internal structure of 
planets and stars, most of the observed astrophysical 
bodies possess a coherent large scale magnetic field. It 
is widely accepted that these natural magnetic fields are 
self-sustained by dynamo action Although reversals 
of the magnetic field in planetary and stellar dynamos are 
now considered to be a common feature, they still remain 
poorly understood. Whereas the Sun shows periodic os- 
cillations of its magnetic field, the polarity of the Earth's 
dipole field reverses randomly. During the last decades, 
several mechanisms have been proposed for geomagnetic 
reversals, among which we can mention, the analogy with 
a bistable oscillator a mean- field dynamo model Q, 
or interaction between dipolar and higher axisymmetric 
components of the magnetic field 0, @] ■ The comprehen- 
sion of dynamo reversals have also benefited from direct 
numerical simulations of the MHD equations, which have 
displayed several possible mechanisms for reversals 0, 

Reversals of a dipolar magnetic field have also been 
reported in the VKS (Von Karman Sodium) dynamo 
experiment In this experiment, periodic or chaotic 
flips of polarity can be observed depending on the mag- 
netic Reynolds number. Based on these results, a model 
for reversals has been recently proposed by Petrelis and 
Fauve ll| . It relies on the interaction between the dipo- 
lar and the quadrupolar magnetic components, and de- 
scribe transitions to periodic oscillations or randomly re- 
versing dynamos. It has been claimed that such a mech- 
anism could apply to the reversals of the Earth magnetic 
field [lH, and temptatively be connected to the periodic 
oscillations of the Solar dynamo. Unfortunately, as for 
many other models, the lack of observations of the mag- 
netic field during a reversal limits a direct comparison 
with the actual geomagnetic reversals. From this point 
of view, the VKS experiment is a unique opportunity to 
test the validity of different models of turbulent revers- 
ing dynamos. In particular, the model makes predic- 
tions about the dynamics that are easily confrontable to 
experimental results. We propose a simple way to ana- 
lyze data from the VKS experiment in order to test this 
model. To wit, we extract from the data the dipolar and 



FIG. 1: (Color online) Sketch of the VKS experiment 



the quadrupolar components of the magnetic field. We 
show that the characteristics of the reversals in the VKS 
experiment are in very good agreement with the predic- 
tions, and that the dynamics of the magnetic field in this 
turbulent dynamo mainly result from an interaction be- 
tween dipolar and quadrupolar modes. 

In the VKS experiment, a turbulent von Karman flow 
of liquid sodium is generated inside a cylinder by two 
counter-rotating impellers, with independent rotation 
frequencies F\ and F 2 (see figure [TJ and [3] for the de- 
scription of the set-up). When Fx = F%, the system is 
invariant about a rotation of an angle n around any axis 
located in the mid-plane. On the contrary, if the im- 
pellers rotate at different rates, this symmetry, hereafter 
referred to as the symmetry, is broken. The dynamics 
of the magnetic field observed in the experiment strongly 
depend on this symmetry. When F\ — F 2 , a statistically 
stationary magnetic field with either polarity is gener- 
ated, with a dominant axial dipolar component. Dynam- 
ical regimes, including periodic oscillations and chaotic 
reversals of the magnetic field, arc observed only when 
the Rk symmetry is broken (F\ ^ F 2 ). 

Previous studies have suggested that the evolution of 
the magnetic field in the VKS experiment results from 
low dimensional dynamics, involving only a few modes in 
interaction ■ This can be ascribed to the proximity of 



the bifurcation threshold and to the smallness of the mag- 
netic Prandtl number in liquid metals (Pm < 10~ 5 ). In- 
deed, in the low- Pm regime, the magnetic field is strongly 
damped compared to the velocity field. Hence, the dy- 
namics are governed only by a small number of magnetic 
modes. Based on this observation, Petrelis and Fauve 



111 ] have proposed that close to the dynamo threshold, 
the magnetic field can be decomposed into two compo- 
nents: 



B = D(t)d(r)+Q(t) q(r) 
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state disappears through a saddlc-node bifurcation [11 
Then, the frequency of the limit cycle vanishes at onset. 
In the vicinity of this transition, turbulent fluctuations 
drive random reversals of the magnetic field. As a 
consequence, random reversals always occur at the 
borderline between stationary and oscillatory dynamos. 
This simple mechanism also yields several predictions 
about the shape of the reversals. First, when the dipolc 
D vanishes, part of the magnetic energy is transferred to 
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where D (respectively Q) represents the amplitude of 
dipolar d(r) (resp. quadrupolar q(r)) component of the 
field. As emphasized in [llj, these components do not 
only involve a dipole or a quadrupole, but also all the 
higher components with the same symmetry in the trans- 
formation R n . In other words, d(r) (resp. q(r)) is the 
antisymmetric (resp. symmetric) part of the magnetic 
field. 

The evolution equations for D(t) and Q(t) can then 
be obtained by symmetry arguments (see [U( for a de- 
tailed description of the model). Since d — > d and 
q — > q in the transformation i?^, dipolar and quadrupo- 
lar modes cannot be linearly coupled when F\ = F2. 
Breaking the R n symmetry by rotating the impellers at 
different speeds allows a linear coupling between dipo- 
lar and quadrupolar modes. For a sufficiently strong 
symmetry-breaking, this coupling can generate a limit 
cycle that involves an energy transfer between dipolar 
and quadrupolar modes. This mechanism has been re- 
cently validated on a numerical model of the VKS exper- 
iment II 511 and also in the case of a mean-field a 2 dynamo 
model [17| . 

Two scenarios of transition from a stationary dy- 
namo to a periodically reversing magnetic field can be 
described in the framework of this dipole-quadrupolc 
model. When the coupling is such that the system is 
close to both a stationary and a Hopf bifurcation, i.e. 
in the vicinity of a codimension-two bifurcation point, 
one can have bistability between a stationary and a 
time periodic reversing dynamo 



We thus get a 
subcritical transition from a stationary dynamo to a 
periodic one with a finite frequency at onset. Turbulent 
fluctuations can generate random transitions between 
these two regimes [l4j |. Far from this codimension-two 
point, a reversing magnetic field can be generated 
through an Andronov bifurcation when the statio nary 
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FIG. 2: (Color online) Parameter space: (cross) no 
dynamo, (circle) stationary dynamos, (star) oscillatory 
or random reversing dynamos. 



the quadrupole Q. An overshoot of the dipolar ampli- 
tude is expected after each reversal. Random reversals 
are asymmetric. During a first phase, fluctuations push 
the system from the stable solution to the unstable one, 
thus acting against the deterministic dynamics. This 
phase is slow compared to the one beyond the unstable 
fixed point, where the system is driven to the opposite 
polarity under the action of the deterministic dynamics. 

In this paper, we use data of the VKS experiment in or- 
der to reconstruct the dipolar and the quadrupolar parts 
of the magnetic field and study their behavior in the time 
dependent regimes. Time-dependent regimes only occur 
for specific values of F\ and F2, inside three delimited 
regions of the parameter space [14j . We will focus on the 
regimes observed when following the line F1/F2 = 0.6 
in the parameter space. Figure [2] shows that when the 
rotation rates are increased along this line, one first bi- 
furcates to a stationary dynamo, then to time-dependent 
regimes. Figure [IJtop) shows the time-recordings of the 
three components of the magnetic field close to the fastest 
disk, displaying the bifurcation from stationary to time- 
dependent dynamo when the frequencies of the two disks 
arc increased from 14.4/24 Hz (Ft + F 2 — 38.4 Hz) to 
15/25 Hz (Fi + F 2 = 40 Hz). After a short transient 
state, the three components of the magnetic field undergo 
a transition to nearly periodic oscillations. 

However, it is hard to test the pertinence of the 
model [ll| from the time recording of the magnetic 
field at a single point. Using measurements obtained 
from two probes 1 and 2, symmetric with respect to 
the mid-plane, we compute the dipolar part, D(i)di = 
(B i (l,t) + B i (2 1 tj)/2 and the quadrupolar part, Q(t)qi = 
(Bi(l,t) — Bi(2,t))/2. In order to obtain observables 
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FIG. 3: (Color online) Top: time recordings of the three 

components of the magnetic field. The rotation 
frequencies are increased from 14.4/24 Hz to 15/25 Hz, 
leading to a transition from a stationary low amplitude 
dynamo regime to a limit cycle. Bottom: Behavior of 
the dipolar and quadrupolar parts of the magnetic field. 
Note the transfer between the two components during 
reversals. 



which arc independent of the spatial component i, each 
of these vectors is projected on its value at a given time 
to- We thus extract D{t) and Q(t) up to a multiplicative 
constant. In the measurements displayed here, the dipo- 
lar and quadrupolar components are projected on their 
stationary values obtained at F\ + F 2 = 38.4 Hz. Note 
however that different methods could be used to recon- 
struct these amplitudes. In particular, plotting the sum 
and the difference of a given component does not change 
the qualitative behavior [lj| ■ FigurcEJbottom) shows the 
evolution of D(t) and Q(t) during periodic oscillations of 
the magnetic field. We observe that when the dipole van- 
ishes, the amplitude of the quadrupole reaches its max- 
imum. This shows that the field reversals observed in 
the VKS experiment do not correspond to a vanishing 
magnetic field, but rather to a change of shape from a 
dominant dipolar field to a quadrupolar one. Immedi- 
ately after each reversal, one can also note that during 
its recovery, the dipolar amplitude strongly overshoots its 



mean value. Therefore, in agreement with the model [ll| . 
reversals in the VKS experiment involve a strong com- 
petition between dipolar and quadrupolar components of 
the magnetic field. 

The decomposition between dipolar and quadrupolar 
components is not only relevant to study these oscilla- 
tions but is also useful to follow the bifurcations observed 
along the line F\/F 2 = 0.6 in figure [2] We now investi- 
gate the evolution of the dynamics in the phase space 
(D, Q) displayed in figured] as Fi + F 2 is modified. The 
limit cycle described in figure [3] bifurcates from a low 
amplitude stationary dynamo when F\ + F 2 is increased 
from 38.4 to 40 Hz. This limit cycle is shown in green in 
figure l4al When F\ + F 2 is decreased again to 38.4 Hz, a 
smaller amplitude limit cycle is obtained (orange curve, 
circles) instead of a fixed point. We need to decrease the 
rotation frequencies further to recover the low amplitude 
stationary dynamo (black dot). Therefore, this transi- 
tion is hysteretic and within some frequency range we 
have bistability involving stationary and time-periodic 
dynamos. This oscillation appears at finite amplitude 
and finite period [3]. The oscillation of figure [3] displays 
a slowing down in the vicinity of two symmetric fixed 
points, as expected for a system close to the saddle-node 
bifurcation of Andronov type. Note however that the on- 
set of the cycle when F\ + F 2 is increased, does not corre- 
spond to such a saddle-node bifurcation, since these two 
stagnation points are distinct from the low amplitude sta- 
tionary dynamo regime obtained at Fl + F2 = 38.4. In 
fact, this transition from a low amplitude stationary mag- 
netic field to an oscillatory regime at finite period, rather 
corresponds to the model taken close to its codimension- 
two bifurcation point [l3|. 

When Fi + F 2 is increased further, the amplitude of 
the limit cycle continuously increases ffigurelibf. In addi- 
tion, the system slows down in the vicinity of two points 
(±D S ,±Q S ) (figure l4c|). Thus, the period of the limit 
cycle significantly increases. For F\ + F 2 = 44 Hz, the 
systems stops on one of these two fixed points and we get 
a stationary dynamo (although we cannot rule out the 
occurrence of other reversals with a longer experiment). 
As explained in the framework of the model [HI], this 
second transition corresponds to a saddle-node bifurca- 
tion or more precisely an Andronov bifurcation: the sta- 
ble fixed points (±D S ,±Q S ) collide with unstable fixed 
points (±D U ,±Q U ) when F\ + F 2 is decreased and dis- 
appear. A limit cycle is thus created, and its period is 
expected to diverge in the vicinity of the saddle-node bi- 
furcation. Turbulent fluctuations of course saturate this 
divergence by kicking the system away from the points 
(±D S ,±Q S ) where it slows down. They also strongly 
modify the dynamics on the other side of the bifurca- 
tion. Indeed, when the stable and unstable fixed points 
are very close one to the other, turbulent fluctuations can 
randomly drive the system from a stable fixed point to its 
neighboring unstable one, and thus trigger a reversal of 
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FIG. 4: (Color online) (a): Evolution of the magnetic field in the phase space (D,Q) at low frequencies: there exists 
a range of bistability in which a stationary dynamo (black dot) and a periodic limit cycle (orange circles) are both 
metastable. (b): Evolution of the limit cycle when F\ + F2 is increased from 40 to 41.8 Hz. (c): Chaotic reversals 
obtained for large values of F\ + F2 in the vicinity of a saddle-node bifurcation. 
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FIG. 5: (Color online) Time evolution of the magnetic 

field during a regime of chaotic reversals, for 
Fi + i*2 = 42.4 Hz. Inset: zoom on the reversal occuring 
at t = 231s. 



the magnetic field. Therefore, random reversals are ex- 
pected in the vicinity of the saddle- node bifurcation [ll[ . 
This is what is observed here as shown below. 

Figure [5] displays the time-recordings of the dipolar 
and quadrupolar components for F\ + F% = 42.4 Hz. 
We observe that both components fluctuate around 
constant values as if they have reached a stable fixed 
point. The time spent in both polarities is random but 
much longer than the magnetic diffusion time scale (of 
order 1 s). One also clearly observes that the amplitude 
of the dipole slowly decreases before rapidly changing 
sign. In the phase space (D,Q) displayed in figure l4cl 
this slow decay corresponds to random motion in the 
regions in the form of elongated spots located along the 
limit cycle in the vicinity of the fixed points. Indeed, the 
motion from each stable fixed point to the neighboring 
unstable one, occurs under the influence of fluctuations 



acting against the deterministic dynamics. It is thus a 
slow random drift compared to the fast reversal phase 
driven by the deterministic dynamics once the system 
has been pushed beyond the unstable fixed point. 
Figure 2c] also show that the spots become more and 
more elongated when F\ + F2 is increased. This tells 
that the distance between each stable fixed point and its 
unstable neighbor increases. Correspondingly, reversals 
are less frequent. For F\ + F2 = 44 Hz (red (light 
grey) cycle in figure I4"c|) , fluctuations can hardly drive 
reversals. As for periodic oscillations obtained above the 
Andronov bifurcation, the modal decomposition (D, Q) 
underlines the short transfer from an axial dipole to 
a quadrupolar magnetic field during random reversals 
obtained below the bifurcation threshold (see inset of 
figure [5]). The dipolar amplitude displays the expected 
behavior, characterized by a slow decay followed by a 
rapid recovery, and showing a typical overshoot after 
each reversal. Evolution in the phase space (D, Q) 
also illustrate how the transfer between dipolar and 
quadrupolar components yields very robust cycles, 
systematically avoiding the origin B = 0. 

In conclusion, we have used a simple method to extract 
the dipolar (antisymmetric) and quadrupolar (symmet- 
ric) components of the magnetic field in the VKS exper- 
iment. We have shown that this decomposition allows 
to investigate the morphology of the magnetic field dur- 
ing reversals, and to compare experimental results to the 
predictions of a recent model proposed in 11] . We have 
shown that the results of the VKS experiment are in very 
good agreement with these predictions: 

- reversals are characterized by a strong transfer to the 
quadrupole when the dipole vanishes, 

- the dipolar mode systematically displays an overshoot 
after each reversal, 

- random reversals are asymmetric, i.e. involve two 
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phases: a slow one triggered by turbulent fluctuations 
followed by a fast one mostly governed by the determin- 
istic dynamics. 

This agreement between the VKS experiment and the 
model has significant consequences. It first shows that 
a fluid dynamo, even generated by a strongly turbu- 
lent flow, can exhibit low dimensional dynamics, involv- 
ing mostly dipolar and quadrupolar modes. Further- 
more, because such a model is based on symmetry ar- 
guments, the mechanisms described here are expected 
to apply beyond the VKS experiment. For instance, al- 
though 3-dimcnsional simulations do not involve a sim- 
ilar level of turbulence, a transfer between dipole and 
quadrupole during reversals has been observed in sev- 
eral numerical studies of the geodynamo @, Q. This 
is consistent with indirect evidences from paleomagnetic 
measurements, suggesting a dipole-quadrupole interac- 
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tion [4J and asymmetric reversals |18j. Observations of 
the Sun's magnetic field also su gge st a transfer between 
dipolar and higher components [16j |. In numerical simu- 
lations based on the VKS experiment [1], a good agree- 
ment with the three predictions reported here has been 
obtained, but only when the magnetic Prandtl number 
is sufficiently small. In this context, our simple method 
could be used to investigate data from numerical simula- 
tions of the geodynamo at low magnetic Prantl number. 
This opens new perspectives to understand the dynamics 
of planetary and stellar magnetic fields with a simple and 
low dimensional description. 

I aknowledgc my colleagues of the VKS team with 
whom the experimental data used here have been ob- 
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